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Abstract. Kasraoui, Stanton and Zeng, and Kim, Stanton and Zeng intro- 
duced certain g-analogues of Laguerre and Charlier polynomials. The moments 
of these orthogonal polynomials have combinatorial models in terms of cross- 
ings in permutations and set partitions. The aim of this article is to prove 
simple formulas for the moments of the g-Laguerre and the g-Charlier polyno- 
mials, in the style of the Touchard-Riordan formula (which gives the moments 
of some q-Hermite polynomials, and also the distribution of crossings in match- 
ings). 

Our method mainly consists in the enumeration of weighted Motzkin paths, 
which are naturally associated with the moments. Some steps are bijective, 
in particular we describe a decomposition of paths which generalises a pre- 
vious construction of Penaud for the case of the Touchard-Riordan formula. 
There are also some non-bijective steps using basic hypergeometric series, and 
continued fractions or, alternatively, functional equations. 



1. Introduction 

Our motivation is to derive in a uniform way generating functions for matchings, 
set partitions and permutations refined by the number of crossings. We achieve 
this by enumerating certain weighted Motzkin paths, which in turn prompt us to 
consider these counts as moments of certain families of orthogonal polynomials. In 
some cases, formulas for these moments are already known. However, the method 
of proof we present in this algorithm is quite general, and leads to very simple 
formulas. 

Let us first define the notion of crossings in matchings, set partitions and per- 
mutations. To do so, it is best to draw the objects we are interested in in a certain 
standard way. We begin with the set of matchings (or fixed-point free involutions) 
A42n of {1, ... , 2n}: these are drawn by putting the numbers from 1 to 2n in this 
order on a straight line, and then connecting paired numbers by an arc. Of course, 
arcs are always drawn in a way such that any two arcs cross at most once, and no 
more than two arcs intersect at any point, see the first picture in Figure [1] for an 
example. Then, a crossing in a matching is, as one would expect, a pair of matched 
points {i,j} and {fc, /} with i < k < j < I, pictorially: 



i k j I 

Indeed, the motivating example for this article is the Touchard-Riordan formula, 
which gives, for each n, the generating polynomial according to crossings for perfect 
matchings of the set {!,..., 2n}. Denoting by cro(M) the number of crossings of 
the matching M, we have: 
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Figure 1. A matching of {1,.. 
associated "histoire de Hermite" . 



10} with 3 crossings and the 
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This has been proved in the 1950's by Touchard, although, curiously, it seems 
that the formula was not given explicitly. This was later rectified by Riordan. 

Quite similar to matchings, a set partition can be depicted by connecting the 
numbers on the line which are in one block B = {bi < 62 < • • ■ < by arcs 
(61,62), (62,63),..., (6/_i,6/), see Figure [2] for an example. Again, a crossing in 
a set partition is what one would expect: a pair of arcs {i,j} and {k,l} with 
i < k < j < I. Denoting the set of set partitions of {1, . . . , n} by n„, the number 
of crossings in a set partition tt by cro(7r) and the number of its blocks by |7r|, we 
will obtain the following g-analogue of the Stirling numbers of the second kind: 



Theorem 1.2. 

(2) 



|7r|=fc 
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where 



[fe] g ~ Y[i=i ^" [i]^^^'' q-binomial coefficient, and [n], 



There is an alternative notion of crossings for set partitions, due to Ehrenborg 
and Readdy [5|, coming from juggling patterns. Namely, we introduce an additional 
infinite arc extending to the right from every maximal element of each block, in- 
cluding singletons, see Figure |3] for an example. Denoting the number of crossings 
in such a drawing of a set partition tt by cro* (tt) , we have: 



Theorem 1.3 (Gould [Qj) 
(3) 
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1 2 3 4 5 6 7 8 




Figure 2. A set partition of {1, . . . , 8} into 4 blocks with 2 cross- 
ings and the associated "histoire de Charlier" . 




1 2 3 4 5 6 7 8 




Figure 3. A set partition of {1, . . . , 8} into 4 blocks with 6 cross- 
ings and the associated "histoire de Charlier-*". 



This is not a new result: essentially, this formula was already known to Gould 
from another definition (the link with crossings is more recent as will appear below) . 

Finally, to depict a permutation cr, we connect the number i with a(i) with an 
arc above the line, if i < <7{i), otherwise with an arc below the line, as done in 
Figure 21 The notion of crossing in a permutation was introduced by Corteel [T] , 
and is slightly less straightforward: a pair of numbers (i, k) constitutes a crossing 
in a permutation^ ii i < k < a{i) < <T{k) or (T(i) < <T{k) < i < k: 

/^?^ or /"VA or \ ^ -^-^--1^ ^ 

i k 7r(i) Tv{k) i k = a{i) a{k) 



Denoting the set of permutations of {1, . . . , n} by ©„, and the number of weak 
exceedances, i.e. numbers i with a(i) > i, of a permutation a by wex(tT), we have: 
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Figure 4. A permutation of {1, . . . , 8} with 5 weak exceedances 
and 5 crossings and the associated "histoire de Laguerre" . 

Theorem 1.4 (Josuat- Verges [12], Corteel, Josuat- Verges, Prellberg, Rubey [2]). 
(4) 



wcx((7) ^cro(cr) 



^ ^' fc=0 \j=0 I \j=0 



yiqi(k+l-i) 



This theorem recently found a rather different proof by the first author '12' . In 
the present article we provide an alternative, using a bijective decomposition of 
weighted Motzkin paths that gives a natural interpretation for the two inner sums. 

The rest of this article is organised as follows. In Section [21 we present some 
background material concerning the combinatorial theory of orthogonal polyno- 
mials. In Section [3J we describe the decomposition of weighted Motzkin paths 
mentioned above, in full generality. Each Motzkin path will be decomposed into a 
Motzkin prefix and another Motzkin path satisfying certain additional conditions. 
In Sectiondl we enumerate Motzkin prefixes, and in Section[5]the other set of paths 
appearing in the decomposition are enumerated. 

There are three appendices. In the first appendix we give an alternative point 
of view of the decomposition presented in Section |3l using inverse relations. In 
the second appendix, we give a bijective proof of the formula for the generating 
function of the paths appearing in the decomposition in the case of set-partitions, 
using a sign-reversing involution. It is thus possible to give a fully bijective proof 
of Theorem II. 2[ analogous to Penaud's proof of the Touchard-Riordan formula. 
Finally, in the last appendix we sketch a proof showing that one cannot expect 
closed forms for Motzkin prefixes with weights different from those considered in 
Section m 



2. Orthogonal Polynomials, moments and histoires 

Motzkin paths are at the heart of the combinatorial theory of orthogonal polyno- 
mials, as developed by Flajolet [5] and Viennot [26;. This theory tells us, that the 
moments of any family of orthogonal polynomials are given by a certain weighted 
count of Motzkin paths. 
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More precisely, by Favard's theorem, any monic sequence of orthogonal polyno- 
mials {Pn)n>o satisfies a three term recurrence of the form 

xPnix) = Pn+l{x) + bnPn{x) + A„P„_i(x), 

where &„ and A„ do not depend on x. Given this recurrence, the n^^ moment fj,^ 
of P can be expressed as the weighted sum of Motzkin paths of length n, that is, 
paths taking up (/^), down (\) and level (— >■) steps, starting and ending at height 
0, and not going below this height, where a horizontal step at height h has weight 
bh and a down step starting at height h has weight A^. 

2.1. Histoires. Three basic examples of families of orthogonal polynomials are 
given by (rescalings of) the Hermite, Charlier and Laguerre polynomials, where the 
moments count matchings (fe„ = 0, A„ = n), set partitions (6,i = 1 + ri, A„ = n) and 
permutations (&„ = 2n + 1, A„ = n^) respectively. It turns out that the Hermite, 
Charlier and Laguerre polynomials indeed have beautiful g-analogues such that the 
moments count the corresponding objects, and q marks the number of crossings. 
We want to establish this correspondence via "histoires" : 

Definition 2.1. Consider a family of orthogonal polynomials with coefficients &„ 
and A„, and fix a„ and c„ such that A„ = a„_ic„ for all n. Suppose that for every 
fixed n, the coefficients a„, 6„ and c„ are polynomials such that each monomial has 
coefficient 1 (as will appear shortly, this is general enough in our context). 

We then call a weighted Motzkin path histoire, when the weight of an up step 
(respectively a level step — 7> or a a down step \) starting at level h is one of the 
monomials appearing in (respectively or Ch)- 

We want to consider four different families of "histoires", corresponding to q- 
analogues of the Hermite, Charlier and Laguerre polynomials. 

Proposition 2.2. There are weight-preserving hijections between 

• matchings M with weight and "histoires de Hermite" defined by 
&„ = 0, a„ = 1 and c„ = [n]q, 

• set partitions tt with weight "histoires de Charlier" defined 
by hn^y+ [n]g, an = y and Cn = [n]q, 

• set partitions tt with weight y^'^^q"° and "histoires de Charlier-^" de- 
fined by bn = yq" + [n]q, a„ = yg" and Cn = [n]q, and 

• permutations a with weight yWox((T)gcro(cr)^ "histoires de Laguerre" de- 
fined by bn = j/g" + [n\q, a„ = yq" and c„ = [n\q. 

These bijections are straightforward modifications of classical bijections used 
by Viennot p6]. We detail them here for convenience, but also because of their 
beauty. . . Examples can be found in Figures HHH 

Proof. The bijection connecting matchings and "histoires de Hermite", such that 
crossings are recorded in the exponent of g, goes as follows: we traverse the match- 
ing, depicted in the standard way, from left to right, while we build up the Motzkin 
path step by step, also from left to right. For every arc connecting i and j with 
i < j, we call i an opener and j a closer. When we have traversed the matching 
up to and including number £, we call the openers i < £ with corresponding closers 
j < £ active. Openers are translated into up steps with weight 1. Accordingly, 
when we encounter a closers £ it becomes a down step with weight q'^, where k is 
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the number of active openers between t and the opener corresponding to It is a 
enjoyable exercise to see that this is indeed a bijection, and that a matching with 
k crossings corresponds to a Dyck path of weight q*' . 

The bijection between set partitions and "histoires de Charher" , due to Anisse 
Kasraoui and Jiang Zeng |14) . is very similar: in addition to openers and closers, 
which are the non-maximal and non-minimal elements of the blocks of the set par- 
tition, we now also have singletons, which are neither openers nor closers. Elements 
that are openers and closers at the same time are called transients. Non-transient 
openers are translated into up steps with weight y, and singletons are translated 
into level steps with weight y. Non-transient closers £ are translated into down 
steps both with weight , where k is the number of active openers between £ and 
the opener corresponding to £. Finally, transient closers I become level-steps with 
weight q^ , with k as before. 

To obtain a "histoire de Charlier-*" of a set partition, using the modified defi- 
nition of crossings, we only have to multiply the weights of steps corresponding to 
closers and singletons by q^ , where k is the number of crossings of the infinite arc 
with other arcs. 

It remains to describe the bijection between permutations and "histoires de La- 
guerre" , due to Dominique Foata and Doron Zeilberger, which is usually done in a 
different way than in what follows, however. To obtain the Motzkin path itself, we 
ignore all the arcs below the line and also the loops corresponding to fixed points. 
What remains can be interpreted as a set partition, and thus determines a Motzkin 
path. Moreover, the weights of the down steps are computed as in the case of set 
partitions, except that the weight of each of those steps needs to be multiplied by y. 
The weights of the level steps that correspond to transients of the set partition are 
also computed as before, but are then multiplied by yq. Level steps that correspond 
to fixed points of the permutation get weight y. The weights of the remaining steps 
are computed by deleting all arcs above the line, and again interpreting what re- 
mains as a set partition. However, this set partition has to be traversed from right 
to left, and weights are accordingly put onto the up steps of the Motzkin path. 
Later, it will be more convenient to move the factor y that appears in the weight 
of all the down steps onto the weight of the corresponding up steps, see Figure [4] 
for an example. □ 

2.2. Particular classes of orthogonal polynomials. In this section we relate 
the families of orthogonal polynomials introduced via their parameters 5„ and A„ in 
Section [2TT] to classical families. We follow the Askey- Wilson scheme [16] for their 
definition. 

The continuous q-Hermite polynomials Hn = Hn(x\q) can be defined |161 Sec- 
tion 3.26] by the recurrence relation 

2x//„ = i/„+i + (l-g")F„_i, 

with H{) = 1. 

Theorem 2.3 (Ismail, Stanton and Viennot [II ). Define rescaled continuous q- 
Hermite polynomials Hn — Hn{x\q) as 

(5) HniAq) = (1 - qr''''Hn{x^\q). 
They satisfy the recurrence relation 

(6) xHn = Hn+l + [n]qHn-l, 
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and their even moments are given by 

(7) E 

The odd moments are all zero. 

The Al-Salam-Chihara polynomials Q„ = (5n(a;; a, can be defined [16[ Sec- 
tion 3.8] by the recurrence relation 

2a;g„ = Qn+i + {a + h)q^Q,, + (1 - - a5g"-i)Q„_i, 

with Qo = 1- We consider two different specialisations of these polynomials. The 
first was introduced by Kim, Stanton and Zeng |15| , and in their Proposition 5 they 
also gave a formula for the moments. However, the formula that follows from our 
Theorem II. 21 appears to be much simpler. 

Theorem 2.4 (Kim, Stanton, Zeng |15j). Define q-Charlier polynomials Cn = 
Cn{x\y\q) as 

(8) c„(.; v\<D = (t^)"'' Qn ( {--y-Tk)-^ 1 1) ■ 

They satisfy the recurrence relation 

(9) xCn = Cn+1 + {y+ [n]q)Cn + y[n]qCn-l 

and their moments are given by 

(10) E 

7ren„ 

The other specialisation was introduced by Kasraoui, Stanton and Zeng |13j . 
however, without providing a formula for the moments (these are actually a par- 
ticular case of octabasic g-Laguerre polynomials from [23]). 

Theorem 2.5 (Kasraoui, Stanton, Zeng |13]). Define q-Laguerre polynomials = 
Ln{x;y\q) as 

(11) Lnix;y\q) = (^)"q„ [(^^2^,^^^^ | . 
They satisfy the recurrence relation: 

xLn = Ln+i + i[n]q + y[n + 1],)Z„ -|- y[n]^L„_i. 
and their moments are given by 

(12) A'n = E 

The Al-Salam-Carlitz I polynomials uli'\x\q) can be defined [El Section 3.24] 
by the recurrence relation 

(13) xUt\x\q) = ul^lMl) + (« + l)9"^l'^H^k) - <r-'a{l - q-)Ul:lM<l), 
with Ul^\x\q) = 1. 
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Theorem 2.6 (de Medicis, Stanton, White [4]). Define modified q-Charlier poly- 
nomials C* = C*{x;y\q) as: 

(14) C:ix; y\q) ^ y^ui^^ (| - ^\q) . 
They satisfy the recurrence relation 

(15) xC: = + (yg" + [n]^)C: + y [n]^ q''-'C:_„ 
and their moments are given by 

(16) /.f = ^ 

7ren„ 

The result from [i] was actuahy stated with another statistic, but both corre- 
spond to Carhtz' q-analogue of the Stirhng numbers of the second kind S[n,k], 
which are such that S[n, k] — S[n — 1. fc — 1] + [fc]gS'[n — 1, fc], and 

n 

/if =^5[n,fc]y^ 

k=l 

3. Penaud's decomposition 

Let us first briefly recall Penaud's strategy to prove the Touchard-Riordan for- 
mula for the moments of the rescaled continuous q-Hermite polynomials Hn. As 
already indicated in the introduction, his starting point was their combinatorial 
interpretation in terms of weighted Dyck paths, down steps starting at level h > 1 
having weight [h]q, up steps having weight 1. 

As the total number of down steps in these paths is n, we may take out a factor 
(1 — <z)^", and instead consider paths with down steps having weight 1 — q^, or, 
equivalently, consider paths with down steps having weight 1 or —q^. 

The next step is to (bijectively) decompose each path into two objects: the first 
is a left factor of an unweighted Dyck path of length n and final height n — 2fc > 0, 
for some k. The second object, in some sense the remainder, is a weighted Dyck 
path of length k with the same possibilities for the weights as in the original path, 
except that peaks (consisting of an up step immediately followed by a down step) 
of weight 1 are not allowed. This decomposition will be generalised in Lemma 13.21 
below. 

The left factors are straightforward to count, the result being the ballot numbers 
(rf-fe) ~ (n-fc-i)" ^'^"^ remainders, Penaud presented a bijective proof that the 
sum of their weights is given by {—l) q^ ^ ' . Summing over all fc we obtain the 
Touchard-Riordan formula ([1} . 

3.1. The general setting. 

Definition 3.1. Let A^„(a, 6, c, d\ q) be the set of weighted Motzkin paths of length 
n, such that the weight of 

• an up step starting at level h is either 1 or — 

• of a level step starting at level h is either d or (a + h)q^ , 

• a down step \ starting at level h is either c or —abq^^^. 
Furthermore, let A^,* (a, 6, c; g) C A^„(a, 6, c, d; g) be the subset of paths that do 
not contain any 

• level step -> of weight d, 
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• peak /^\ such that the up step has weight 1 and the down step has weight 
c. 

Finally, let 'Pn,k{c,d) be the set of left factors of Motzkin paths of length n and 
final height k, such that the weight of 

• an up step is I, 

• a level step — >■ is d, 

• a down step \ is c. 

With these definitions, the decomposition used by Penaud can be generalised in 
a natural way as follows; 

Lemma 3.2. There is a bijection A between A^„(a, b, c, d; q) and the disjoint union 
of the sets Vn,k{c, d) x A^^(a, b, c; q) for k G {0, . . . ,n}. 

Proof. Let H be a path in Mn{a, b, c, d; q). Consider the maximal factors fi, . . . ,fj 
of H that are Motzkin paths and have up steps of weight 1 , level steps of weight d 
and down steps of weight c. We can thus factorise H as /10/1/11/2 ■ • ■ fjhj. 

Since this factorisation is uniquely determined, we can define A(_ff) = (Hi,H2) 
as follows: 

i^i = (^)l''°l/i(^)l''^l/2.../,(^)l''^l and H2^ho...hj. 
Thus, Hi is obtained from H by replacing each step in the hi by an up step 
and H2 is obtained from H by deleting the factors f.^. Since the /j are Motzkin 
paths, the weight of H is just the product of the weights of Hi and H2 ■ Furthermore, 
it is clear that Hi is a path in Vn.k{c, d) with final height k = \ho\ + \hi \ + ■ ■ ■ + 
\hj\. We observe that the hi cannot contain a level step — > of weight d or a peak 
such that the up step has weight 1 and the down step has weight c, because 
then the factorisation of H would not have been complete. Thus H2 is a path in 
Ml{a,b, c; q). 

It remains to verify that A is indeed a bijection. To do so, we describe the inverse 
map: let {Hi,H2) € Vn,k{c, d) x A^^(a, 6, c; q) for some fc G {0, . . . , n}. Thus, there 
exists a unique factorisation 

ifi = (^r/i(^)"V2.../,(^r 

such that the fi are Motzkin paths and k = J2i=o Write H2 as /lo • ■ • hj, where 
the factor hg has length ui. Then A~^{Hi, H2) = /io/i^i./2 ■ • ■ fjhj is the preimage 
oi{Hi,H2). ' □ 

3.2. Specialising to matchings, set partitions and permutations. As re- 
marked in the introduction of this section, we begin by multiplying the weighted 
sum of all Motzkin paths by an appropriate power of 1 — g. In the case of matchings 
of {1, ... , 2n}, we are in fact considering Dyck paths of length 2n where a down step 
starting at height h has weight [h]q. Multiplying the weighted sum with (1 — g)", 
or, equivalently, multiplying the weight of each down step by 1 — g, we thus obtain 
Dyck paths having down steps starting at height h weighted by 1 — g'*, which fits 
well into the model introduced in Definition 13.11 namely, the set A^„(a, 6, c, d; 
with a = 0, 5 = 0, c=l and d = consists precisely of these paths - except that 
they are all reversed. 

In the case of set partitions of {1, . . . , n}, multiplying the weighted sum by (1 — 
5)" and reversing all paths we see that we need to enumerate the set A^„(a, b, c, d; q) 
with a = b — —1, c — y{l — q) and d = 1 + y{l — q). When using the modified 
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definition of crossings in set partitions, we obtain surprisingly different parameters, 
namely a = — 1, b = y{\ — q), c = and d = 1. Finally, the case of permutations of 
{1, . . . , n} is covered by enumerating the set A^„(a, 6, 1, d; q) with a = —1, b ~ —yq, 
c = y and d = 1 + y. 

4. Counting Vn,kic,d) 

In general, formulas for the cardinality of 'Pn,k{c,d) can be found easily using 
Lagrange inversion [53]. Consider the generating function Pk = \T^n,k{c, i", 
we want to determine the coefficient of in tPk = {tPo)'^^^ . Observing the 
relationship 

tPo = t {1 + d{tPo) + c{tPof) 



we find that [i"](tPo) = ^[z''-'']{l + dz + 



and thus 



(17) 



\Vn,k{c, d)\ 



71— k 

1 / n 



1 ^ 

1=0 



I 



21- 



^2l — n+k^n—k—[ 



To count matchings, set partitions or permutations according to crossings (mod- 
ified or not), the only sets of parameters that we need to consider are (c, d) — (1, 0), 
(c, d) = (1, 2) and (c, d) = (0, 1). Curiously, these are precisely the values for which 
Equation ()17|) allows a closed form, i.e. can be written as a linear combination of 
hypergeometric terms. A (sketch of a) justification of this fact is given in Appen- 
dix O 

4.1. Matchings. For matchings, we have {c,d) = (1,0) and we obtain the ballot 
numbers; 

Lemma 4.1. The cardinality of Vn.n-2k{^,0) , i-e. the number of left factors of 
Dyck paths of length n and final height n — 2fc > is 



n 
fc- 1 



4.2. Set partitions and permutations. For set partitions and permutations, we 
have (c, d) = {y,l + y) and obtain the following: 

Lemma 4.2. The generating function for Vn.k{y, 1 + J/) is: 

n-k 

(18) E 



n 



n 



n 

j + k- 



1 



y 



Proof. The elements of Vn,k{y, 1 + y) have weight 1 + y on each level step. However, 
it is again more convenient to pretend that there are two different kinds of level 
steps, with weight 1 and y respectively. Let P be a left factor of a Motzkin path 
with weight y^ . We then use the following step by step translation to transform it 
into a pair {Ci, C2) of non-intersecting paths taking north and east steps, starting 
at (0,1) and (1,0) respectively (see Figure [S] for an example): 



i'h step of P 


step of Ci 


i"^ step of C2 


weight 1 


t 
t 




— >, weight y 






\, weight y 




t 
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/ 


1 


y 


y 


1 


/ 







n - j + 1 

n — k ~ j 
1 




1 j j+k+1 
Figure 5. A bijection to count 'Pn,k{y, 1 + y)- 

The condition that a Motzkin path does not go below the a;-axis translates into 
the fact that Ci and C2 are non-intersecting. Since P has j steps weighted by y, 
the path Ci ends at (j, n — j -I- 1). Since P ends at height k, the number of up steps 

and the number of level steps — >■ with weight y add up to j + fc, so C2 ends at 
{j + k + l,n- j - k). 

By the Lindstrom-Gessel-Viennot Lemma [8], these pairs of non-intersecting 
paths can be counted by a 2 x 2-determinant, which gives precisely Formula ([T8| . □ 



For y — 1 the sum in Equation (1181) can be simplified using Vandermonde's 
identity. Thus, the number of left factors of Motzkin paths of length n and final 
height k, with weight 2 on every level step is 



(19) 



2n 

n ^ k 



2n 

-k-2 



For {c,d) — (0, 1), that is, y — 0, we obtain what we need to count modified 
crossings in set partitions, namely the binomial coefficient (^) . 



5. COUNTING Ml{a,b,c;q) 

In this section we use a continued fraction to find the generating function for 
the Motzkin paths in A^^(a, b, c; q) (these paths are described in Definition 13. ip . It 
turns out that this continued fraction can be expressed as a basic hypergeometric 
series, which allows us to compute the coefficients corresponding to paths with given 
length. Let K{a, b, c; q) be 



(20) 



1 



l + c-{a + b) 



ic-ab)il-q) 



(a + b)q 



{c-abq){l-q^) 



1 ^ . ^ M 2 {c-abq^){l-q^) 
1 + c — (a + b)q'^ 



Let us first give a combinatorial interpretation of K{at,bt,ct^]q) in terms of 
weighted Motzkin paths. This result is close to those given by Roblet and Vicn- 
not [22], who developed a combinatorial theory of T- fractions. These are continued 
fractions of the form 1/(1 — aot — bot/ (1 — ait — bit/ ...)), and they are generating 
functions of Dyck paths with some weights on the peaks. 

Proposition 5.1. The coefficient of t^ in the expansion of K{at,bt,ct^;q) is the 
generating function of M'^{a, b, c; q). 
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Proof. The continued fraction K{at, bt, ct^; q) equals: 
(21) 

1 



l + cf^ - {a + b)t- 



t^{c-ab){l - q) 



1 + - {a + b)qt 



t\c~abq){l - q^) 



t^(c-abq^){l-q^) 
l + ct^-{a + b)qH ^ — — 



Using the ideas introduced by Flajolet O Theorem 1], we thus obtain paths with 
four types of steps, denoted up down \, level — )• and double-level — >, the last 
type of step simply being twice as long as the usual level step. Moreover, the weight 
of 

• an up step starting at height h is either 1 or — g''^^, 

• a level step — ^ starting at height /i is (a + b)q*^ , 

• a down step \ starting at height h is either c or —abq^^^, 

• a double-level step — >■ is — c. 

To prove the statement, it suffices to construct a involution on the paths, such that 

• its fixed points are precisely the elements of A^^(a, 6, c; q), i.e. paths with- 
out double-level steps — > and without peaks /^\ such that the up step 
has weight 1 and the down step has weight c, 

• the weight of a path that is not fixed under the involution and the weight 
of its image add to zero. 

Such an involution is easy to find: a path that is not in A^^(a, 6, c; q), we look for 
the first occurrence of one of the two forbidden patterns, i.e. a double level step 
— > or a peak /^\ with steps weighted 1 and c respectively. We then exchange 
one of the patterns for the other - since the double level step — > has weight — c, 
the weights of the two paths add up to zero. □ 

As mentioned above, K{at, bt, ct^; q) can be expressed as a basic hypergeometric 
series. We use the usual notation for these series, as for example in [7]. 



Proposition 5.2. For A ^ 1, B ^ 0, we have 

q,B 

1 — \ 

For A ^ 1, B — Q, we have 

q,Cq 



(22) KiA,B,C;q) = ^ • ^q^i (^^'^''^ 



(23) K{A,0,C;q) = ^^-,^,(^j^ 

Proof. Consider the following more general continued fraction, containing a new 
variable z: 

l + C-iA + B)z iC-ABqz)il-q^z) 

l + C-(A + B)qz ^ ^— Vtt ^ 

. . o iC - ABq^zMl-q^z) 
1 + C - A + B)q^z - ^—^ — 
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Following Ismail and Libis [10] (see also Identity 19.2.11a in the Handbook of Con- 
tinued Fractions for Special Functions [3j), we have: 

-1 



M{z) = 



1 



1 



2<Pl 



A,B 
Cq 



q,qz 



2<Pl 



A,B 
Cq 



q,z 



To be able to specialise z = 1, we can use one of Heine's transformations fP, p. 13]. 
For i? ^ we obtain 



^ 1 {Aqz,B,Cq,z;q)oc 

M{Z) = — ;rr-7:; ^ • 2<Pl 



1 — z (Az, B,Cq,qz;q) 
1 



3 q,qz 
Aqz 



q,B 



fCB~^q,z 
[ Az 



q,B 



fCB-^q,qz 



1- Az 
In case i? = 0, we have 

1 {Aqz,q,Cq,z;q) 



q,B 



291 



^CB -q,z 



Az 



q,B 



M{z) = 



1 — z (Az,q,Cq,qz;q) 
1 



\Aqz 



q,Cq] ■ 10 



z 

Az 



q,Cq 



1 



\Aqz 



q,Cq 



■ l(Pl 



q,Cq 



□ 



Remark. Although the symmetry in A and B is not apparent in Equation (|22l) . it 
can be seen using one of Heine's transformations [71 p. 13]. 



In the following, we will always use 



(24) 



K{at,bt,ct^;q) 



1 



1 - at 



cb ^qt, q 
aqt 



q,bt 



Besides, it is also possible to use a method giving M{z) as a quotient of basic 
hypergeometric series without knowing a priori which identity to use. This method 
was employed in [5], following Brak and Prellberg [T^]. Namely, note that the 
continued fraction expansion of M(z) is equivalent to the equation: 

(25) M(z) = ^ ^ — . 

^ ^ l-c+{a + b)z- {c- abz){l-qz)M{qz) 

By looking for solutions of the form M{z) = (1 — az)~^ ^h(z) ■• obtain a linear 
equation in H{z), which gives a recurrence for the coefficients of the Taylor expan- 
sion oi H{z), which is readily transformed into the explicit form oi H{z) as a basic 
hypergeometric series. 

5.1. Matchings. For matchings, we have (a, &, c) = (0,0, 1) and obtain: 
Lemma 5.3. 

oo 

X(0,0,t2;g) = 5:(-<ygm. 

k=0 

Essentially, this was shown by Penaud 17^, who enumerated A^2fe(0, 0, 1; g) by 
first constructing a bijection with parallelogram polyominoes, passing through sev- 
eral intermediate objects with beautiful names like 'cherry trees'. On the polyomi- 
noes he was finally able to construct a weight-preserving, sign-reversing involution. 



14 MATTHIEU JOSUAT-VERGES AND MARTIN RUBEY 

with the only fixed point having weight (— Ij gv ^ corresponding to weighted 
Dyck paths with a single peak, and all weights maximal. 

5.2. Set partitions. For set partitions, we have (a, 6, c) = (O, — — q)) and 
can use the following lemma: 

Lemma 5.4. 

(26) 

Proof. Using Equation (f24| we find: 

K{0,-t,ct^;q)^ 2(1)1 



1=0 j=0 



-cqt,q 




The proof follows by plugging in the elementary expansion 



{-cqt-q),^\{{l + q^ct) = Y,qi'^^") 



3=0 



C>tK 



□ 



In the appendix, we give a bijective proof of this lemma. 



5.3. Set partitions, modified crossings. When using the modified definition of 
crossings in set partitions, we have (a, 6, c) = ( — 1,2/(1 — ?),0) and can use the 
following lemma: 



Lemma 5.5. 

(27) 

Proof. Using Equation (l24l) we find: 

1 



i^(-t,6t,0;g) =^^f6^(-i: 

1=0 j=0 



K{-tM,Q;q) 



l + t 



2<Pl 



0,q 

-qt 



\ 1 



The proof follows by plugging in the elementary expansion 

1 A 1 



(- 



. i . 



□ 



5.4. Permutations. In the case of permutations, we have (a, 6, c) — {—l,—yq,y) 
and find: 



Lemma 5.6. 



K{-t, -yqt, yt'; q) - £(-i)'= ^ z/g'C'+i- 



A;=0 
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Proof. Using Equation (f24l) . we have: 

" {-yqty 



K{-t,-yqt,ye-q)- ^ W 9 



l + t ~qt 



OO 



q.-yqt] Y 



i=0 

oo oo oo oo 



1=0 j"=0 i=0 j=0 

To finish the proof it only remains to substitute k for i + j. □ 

6. Conclusion 

Let us briefly summarize how the four theorems announced in the introduction 
can be proved using the previous sections. In each case, the enumeration of crossings 
in combinatorial objects is linked with the enumeration of the weighted Motzkin 
paths in Mn{a,b,c,d;q). The bijection A shows that the generating function of 
crossings can be decomposed into the generating functions of the sets Vn,k{c, d) and 
A^* (a, 6, c; g), which in turn have been obtained in the previous two sections. This 
fulfills our initial objective as stated in the introduction. 

Appendix A. Inverse relations 

We would like to mention an interesting non-bijective point of view of the path 
decomposition given in Section|31 using inverse relations. Given two sequences {a„} 
and an inverse relation is an equivalence such as, for example: 

Vn>0, a« = E(^)6fe ^ Vn > 0, = E (^) (-ir^'afc- 

fc=0 ^ ^ k=0 ^ ^ 

This particular relation is easily proved by checking that the (semi-infinite) lower 
triangular matrix ((/))i.jgN has an inverse, which is ((-l)'+J'(*.)),,jeN. Other rela- 
tions of this kind can be found in Chapters 2 and 3 of Riordan's book 'Combinatorial 
Identities' [2]. To prove the Touchard-Riordan formula ([T]), let a,2n — (1 — q)"f^§n 
and a2n+i — 0. We then use the following inverse relation ;21, Chapter 2, Equa- 
tion (12)]: 

LfJ LfJ 

(28) -n-Y.{il)-L-l))bn-2k ^ 6„=E(-l)'(V)«-2fe. 
fe=0 fc=0 

Again, this can be proved by inverting a lower-triangular matrix. It remains to 

prove that b2k = (— I)''?'- ^ ^ and b2k+i — 0. To this end, we relate the sequence 
bn to Schroder paths: 

Definition A.l. A Schroder path of length 2n is a path in N x N starting at (0, 0), 
arriving at (2n, 0) with steps (1, 1), (1, —1), or (2, 0). 

Lemma A. 2. Suppose that 02,1+1 = and a2n is the generating function of Dyck 
paths of length 2n, with weight 1 — on each north-east step starting at height 

h (this is to say a2n = (1 ^ 'Z)'V2riy'- Suppose that a„ and 6„ are related by (j28p . 
Then 62n+i = 0, and b2n is the generating function of Schroder paths of length 2n, 
with weight —1 on each level step, and 1 — q^^"^ on each north-east step starting at 
height h. 
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Proof. For any even n, consider a Schroder path of length n with k level steps, 
weighted as described above. This path has n — k non-level steps, and can thus be 
obtained from a Dyck path of length n — 2fc by inserting the level steps. There are 
("fe '^ays to do so, which implies that the generating function of Schroder paths 

is indeed equal to X)i=o(~-'^)'^("fe'^)'^"-2fe, and therefore equal to 6„. □ 

With the bijective decomposition of paths in Section[31 we showed that we have to 
obtain the generating function of the set A^^.(0, 0, 1; q) to prove Touchard-Riordan 
formula. With the inverse relations, we showed that we have to count certain 
weighted Schroder paths. The fact that both sets have the same generating function 
follows from the involution given in the proof of Proposition 15.11 

It is also possible to use an inverse relation to obtain the formula in Theorem ll.4l 
for the g-Laguerre moments. By inverting a lower triangular matrix, one can check 
that 

(29) - E ( E ((") J - o-i) 0+iVi)) I 

k=0 \j=Q J 

for all n is equivalent to 

^^i:(|\„!;^.)(""rO'-»'''--"H°' 

for all n. 

Equation ((30)) can be interpreted as follows: given that ak counts elements of 
A4n{l, —yq, 1 + y, q), then bn count paths of length n, with the same weights as 
in A^„(l, —yq, y,l + y; q), where we insert some level steps — )• with weight —I — y, 
and some double level steps — > with weight —y. 

Indeed, suppose that we inserted j double- level steps — >, and hence n — fc — 2j 
level steps — starting with a path of length k. This yields the weight (— j/)"' (— 1 — 
y^n~k-2j ^j^g inserted steps. The first binomial coefficient, is the 

number of ways to insert the n — k — j level steps among the n — j steps (the total 
number of steps being n — j because the length is n, and j steps of double length). 

The second binomial coefficient, ("~*^~"'), is the number of combinations of the 
j inserted double-level steps and the n — k — 2j inserted level steps. Using the 
involution given in the proof of Proposition 15.11 we see that 6„ counts elements in 
the set M*^{l,''yq,y;q). 

Appendix B. A bijective proof of Lemma [5T4l 

We show in this appendix that Penaud's bijective method of proving Lemma l5.3l 
can be generalised to prove also Lemma l5.4l Namely, we construct a sign- reversing 
involution on a set of weighted Motzkin paths, whose fixed points are enumerated 
by the right-hand side of (pS]) . This involution was essentially given by the first 
author in [12] in a different context. We take the opportunity to correct some 
mistakes in this reference. 

In the following we fix integers j,k > and consider the set Cj^k of Motzkin 
paths of length k + j with k — j level steps — >, (and hence j up steps and j down 
steps \), satisfying the following conditions: 

• the weight of all up steps is 1, 
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• the weight of a level step — at height h is —q^, 

• the weight of a down step \ starting at height h is either 1 or — g'', 

• there is no peak /'Xj such that both the up step and the down step have 
weight 1. 

The generating function of Cj,k is then the coefficient of a^t^^^ in K{Q, —t, at^; q). 
It thus suffices to prove the following: 

Proposition B.l. There is an involution 9 on the set Cj^k such that: 

• the fixed points are the paths that 

— start with j up steps 

— and contain no down steps \, of weight 1, 

• the weight of a path that is not fixed under the involution and the weight of 
its image add to zero. 

Moreover, the sum of weights of fixed points of the involution is (— l)'"^'- 2 ) , 

J q 

Penaud's method consists in introducing several intermediate objects as de- 
scribed in Section 15.11 However, in the case at hand we will not use intermediate 
objects, but rather construct the involution directly on the paths. What we give is 
a generalisation of Penaud's construction, which we recover in the case k = j . 

Proof. Following Penaud |17) , we use in this proof a word notation for elements in 
Cj^k- The letters x, z, y, and y will respectively denote the steps — \ with 
weight 1, and \, with weight — q''. For any word c G Cj^fc, we define: 

• u{c) as the length of the last sequence of consecutive x, 

• v{c) as the starting height of the last step y, if c contains a y and there is 
no X after the last y, and j otherwise. 

See Figure [S] for an example. The fixed points of 6 will be c G Cj^k such that 
u{c) — v[c) — j, which correspond to the paths described in Proposition lB.il 

now, suppose c is such that u{c) < j or v{c) < j. We will build 6 so that 
v{c) < u{c) if and only if u{9{c)) < v{9{c)). Thus it suffices to define 9{c) in the 
case v{c) < u{c), and to check that we have indeed u{9{c)) < v{9{c)). So let us 
suppose v(c) < u{c), hence v{c) < j. 

Since w(c) < j, there is at least a letter y in c having no x to its right. Let c 
be the word obtained from c by replacing the last y with a y. There is a unique 
factorisation 

c = /ix"(^) 02/^2 

such that: 

• a is either z, or y, 

• /2 begins with z or y, contains at least one letter y, but contains no x. 
Let us explain this factorisation. By definition of u(c), we can write c = /ix^^'^^c', 
where c' does not contain any x. In a word c G Cj^k, an x cannot be followed by a y. 
So we can write c = /i where a is either z, or y. Then, we write c" = 
with £ > maximal, and /2 satisfy the conditions (/2 contains indeed a y because 
we transformed a y into a y). Uniqueness is immediate. 

We set: 

(31) 9ic) = /ix"(^)-''(^)ay^a;''(^)/2. 
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See Figure [6] for an example with u(c) = 4, v{c) = 2, j = 9 and k = 12. We can 
check that w{c) = ~q^^ = ~w{e{c)), and u{e{c)) = 2, u(6'(c)) = 3. 



u{c) = 4 




xxyxzyxxyxxxxzyyyyyzy 



c = 




xxyxzyxxyxxxxzyyyyyzy 




xxyxzyxxyxxxxzyyyyyzy 



Figure 6. An element c in Cj^k, and its image 0{c). Thick steps 
are the steps y and z, i.e. steps with weight —q^. 

We show the following points: 

• The path 6{c) is a Motzkin path. Indeed, the factor ay^ in c ends at height 
at least w(c), since the factor /2 contains a step y starting at this height 
and contains no x. We can thus shift this factor c so that the result is again 
a Motzkin path. 

• The path c and its image 6(c) have opposite weights. To begin, between 
c and c, the weight is multiplied by —q'"'^'^\ since we have transformed a 
y into a y starting at height v{c). Between c and 0(c), the height of the 
factor ay^ has decreased by w(c), so the weight has been divided by q'"''^\ 
A factor —1 remains, which proves the claim. 

• The path 9{c) is such that u{9{c)) < v{9{c)). From the definition (|5T|) we 
see that u{9{c)) — v{c). Besides, v{c) < v{9{c)) since the last step y of c 
has been transformed into a y to obtain c and 9(c). 

• Every path c' with u(c') < v(c') is obtained as a 9(c) for some other path 
c with u(c) > v(c). Indeed, let c' be the word obtained from c' after 
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replacing the last y at height u(c') with a y. There is a unique factorisation 
c' = fiay^x'"^'^ '/2, where a is z or y, and /2 contains no x. Then by 
construction, c = fix'^^'^ ^ ay^ has the required properties. 

Thus, 6 is indeed an involution with the announced fixed points. 

It remains only to check that the sum of weights of the fixed points is equal to 
(— 1)'^(7(^2 ) l^'^j , A fixed point of 9 is specified by the heights /ii, . . . , hk~j of the 
k — j steps — >. These heights can be any set of values provided that j > hi > ■ ■ ■ > 
hk-j > 0. The weight of such a fixed point is: 

fe-i 

Indeed, the j steps have respective weights — g, — g^, . . . , , which gives a factor 
[j-Xyq . Besides, we have: 

E i^'^ = \''' , 

j>hi>--->hk-j>o LJJ, 
by elementary property of g-binomial coefficients. This ends the proof. □ 

Appendix C. Closed forms for \Vn,kic,d)\ 

In this appendix we give a justification of the fact that there is no (hypergeo- 
metric) closed form (in the sense of Petkovsck, Wilf and Zeilberger [18]) for 

(32) |P„,*, .)l ^ ^ g (" I ') _ 1 , . 

except when (c, d) is one of (1,0), (0, 1) or {y^,2y). More precisely, we claim that 
iVn.kic, d)\ cannot be written as a linear combination (of a fixed finite number) of 
hypergeometric terms except in the specified cases. In the following we sketch a 
straightforward way to check this is using computer algebra. 

First we convert the summation into a polynomial recurrence equation. This can 
be done by using Zeilberger's algorithm (which also proves that the recurrence is 
correct), for example. Writing p„ = \Vn,k{c,d)\ we obtain 

(33) {4:C-d^){n+l){n + 2)pn+d{n+2){2n+5)pn+i-{n + 2-k){n+4:+k)p„+2 = 0. 
Alternatively, one can also find a recurrence for qk = iVn.kic, d)\, which is 

(34) (fc + 3)(fc - n)qk + d{k + l)(fc + 3)qk+i + c{k + l){k + n + 4)qk+2 = 0. 

It now remains to show that both equations admit no hypergeometric solutions, 
except for the values of (c, d) mentioned above. To this end we use Petkovsek's 
algorithm hyper, as described in Chapter 8 of 'A=B' [TH]. Unfortunately, this 
time we cannot use the implementation naively. Namely, a priori hyper decides 
only for fixed parameters (c, d) whether a hypergeometric solutions exists or not. 

However, it is possible to trace the algorithm, and, whenever it has to decide 
whether a quantity containing c or c? is zero or not, do it for the computer. (Of 
course, it should be possible to actually program this, but that is outside the scope 
of this article.) We refrain from giving a complete proof, but rather give only a few 
details to make checking easier. 
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First of all, let us assume that c, d and (P ~ 4c are all nonzero. Then the degrees 
of the coefficient polynomials in both recurrence equations are all the same. From 
now on, the procedure is the same for both recurrence equation, so let us focus 
on the one for p„. According to the remark in Example 8.4.2 in |18| . we have to 
consider all monic factors a{n) of the coefficient of p„, and also the monic factors 
b(n) of the coefficient oi pn+2, such that the degree of a(n) and b{n) coincide. In 
this case, the characteristic equation. Equation (8.4.5) in one has to solve turns 
out to be — 2dz + rf^ — 4c. For each of the two solutions in z, one has to check 
that there is no polynomial solution of the recurrence 

Pain)cn + zPi{n)Cn+i + z'^P2{n)Cn+2 = 0, 

where the coefficient polynomials ^0(^^)1 and ^2('t-) sue polynomials derived 

from the coefficient polynomials of the original recurrence by multiplying with cer- 
tain shifts of a{n) and b{n). 

This can be done with the algorithm poly, described in Section 8.3 of |18) . 
Namely, depending on the degrees of yet another set of polynomials derived from 
Po{n),Pi{n) and P2{n), it computes an upper bound for the degree of a possible 
polynomial solution. Indeed, the algorithm decides that the degree of such a solu- 
tion would have to be negative, provided that c, d and d'^ — 4c are nonzero, which 
is what we assumed. 
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